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Abstract 



We study both the branching ratio for 6 — > 37 decay and the muon anomalous 
magnetic moment, a^ = {g — 2)^/2, in the minimal supersymmetric standard 
model with gauge-mediated supersymmetry breaking. Combining new exper- 
imental data on and the branching ratio for b ^ s^, strong limits on the 
parameter space of these models are derived. We find that this combined 
study leads to much stronger constraints on the parameter space of the model 
than those from either 6 — > 57 or a^. In particular, the region of large tan P is 
extremely limited, which would have been otherwise allowed. We include the 
supersymmetric one-loop correction to the mass of b quark, rrn,, and find that 
in order to have a correct value of rub, the region of large tan/3 and /i < (in 
our convention) is not allowed in these models. The region of large tan/3 and 
H > is also strongly constrained. We present bounds on supersymmetric 
particle masses as a function of tan /3. 
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Supersymmetry (SUSY) is a very attractive candidate beyond the standard model (SM) 
since it provides an elegant solution to the hierarchy problem in particle physics. Under- 
standing the mechanism of SUSY breaking and its communication to the observable sector 
is one of the most important open questions. There are two types of realistic supersymmet- 
ric models of interest : gravity-mediated models and models with gauge-mediated SUSY 
breaking (GMSB). The GMSB models have been of special interest, because they have 
attractive features of natural suppression of the SUSY contributions to flavor-changing neu- 
tral currents at low energies and prediction of the supersymmetric particle mass spectrum 
in terms of few parameters. 

The parameters of SUSY models can be constrained by using precision measurements 
in low energy experiments, because superparticles contribute to low energy physics through 
radiative corrections. In particular, experimentally observed rare decays may shed light 
on the parameter space of SUSY models. Processes such as 6 — >• S7 do not occur at the 
tree level, and at one-loop level they occur at a small rate but enough to be sensitive to 
new physics effects 0. The CLEO collaboration has recently reported the branching ratio 
for the decay 6 ^ 57 § : BR(6 ^ 57) = (3.15 ± 0.35 ± 0.32 ± 0.26) x 10-^ which 
corresponds to the bound 2.0 x 10"^ < BR(6 57) < 4.5 x 10"^ at 95 % C.L. The 
anomalous magnetic moment of muon, = {g — 2)^/2, is also sensitive to new physics 
effects and can be used to constrain SUSY models on account of the great accuracy 

of both experimental and SM theoretical values of a^. The present experimental value of 
is a^^P = 11659230(84) x lO"^'', while the theoretical prediction for in the context of 
the SM [1,1 is af^ = 11659162(6.5) x 10"^°. 

In this paper, we obtain combined constraints due to both b ^ s'j decay and in the 
minimal supersymmetric SM (MSSM) with GMSB. Even though there exist the previous 
works which studied either 6 — 57 [|T3,|n[l '^m ^^e GMSB models, our work extends 
the previous ones in the sense that we investigate both b s'j and together with the 
inclusion of the supersymmetric one- loop correction to the mass of b quark, rub, which has 
considerable effects in large tan j3 region |[T2|J13[| . We shall see that this combined study leads 
to much stronger constraints on the parameter space of the model than those from either 
b —>■ s'-f or a^. In particular, the region of large tan/5 is extremely limited, which would 
have been otherwise allowed. Furthermore, in this work, we explicitly show that, with the 
presently available experimental data, constraints from the decay 6 — 57 are more stringent 
than those from in broad region of the parameter space. We also present bounds on 
supersymmetric particle masses as a function of tan /5. 

In the GMSB models messenger fields transmit SUSY breaking to the fields of visible 
sector via loop diagrams involving SU(3)cxSU(2)ixU(l)y gauge interactions. The simplest 
model consists of messenger fields which transform as a single flavor of vectorlike 5 + 5 of 
SU(5). These messenger fields may be coupled to a SM singlet chiral superfield 5* through 
the superpotential 

^messenger = )^dSDD + Xj^SLL, (l) 

where the fields have the SM representations and quantum numbers D : (3, l)y=_2/3, D '■ 
(3, l)y=2/3, L : (l,2)y=_i, and L : (l,2)y=i. The scalar and F components of S acquire 
VEVs (S) and (Fs), respectively, through their interactions with the fields of hidden sector, 
which results in breakdown of SUSY. It is known that for messenger fields in complete SU(5) 
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representation, at most four (5 + 5) pairs, or one (5 + 5) and one (10 + 10) pair are allowed to 
ensure that the gauge couplings remain perturbative up to the grand unified theory (GUT) 
scale H . 



In general, the parameters /i and B in soft SUSY breaking terms depend on the details 
of the SUSY breaking in the hidden sector. We require that electroweak symmetry be 
radiatively broken, which determines /i^ and B in terms of other parameters of the theory. 
Then the sparticle masses depend on five independent parameters: M, A, n, tan/3, and 
sign(/x). Here M is the messenger scale given by M = A (S) with a universal Yukawa coupling 
A in the messenger sector at GUT scale. The parameter A is defined by A = (Fs) / (S). The 
integer number n is the effective number of messenger fields given hj n = + 3nio, where 
ns and nio denote the number of (5 + 5) and (10 + 10) pairs, respectively. The radiatively 
generated soft SUSY-breaking masses of gaugino and scalars, Mj and m^, at messenger scale 
M are given by ||15|,|16 



m{M) = ng{x)^^ A, (2) 



a 



An J ' 



(3) 



where x = A/M. (z = 1, 2, 3) are the three SM gauge couplings with GUT normalization 
for ai. ki are 1, 1, 3/5 for SU(3)c, SU(2)i and U(l)y, respectively. Ci are zero for gauge 
singlets, and 4/3, 3/4 and {Y/2y for the fundamental representations of SU(3)c, SU(2)i, 
and U(l)y, respectively (with Y defined by Q = I3 + Y/2). g{x) and /(x) are messenger 
scale threshold functions. 

We use the input values asiMz) = 0.118, sin^ OwiMz) = 0.2315 and a{Mz) = 1/128. 
The parameter A is taken to be around 100 TeV to ensure that the sparticle masses are of 
the order of the weak scale. We restrict 1 < M/A < 10^ : the case M = A is excluded 



since it produces a massless scalar in the messenger sector [|T5|,|16| and the upper bound on 



the gravitino mass of about 10^ eV restricts M/A < 10^ |[16|J17[| . Using the appropriate 



renormalization group equations (RGEs) |]T8[, we first go up to the messenger scale M with 
gauge and Yukawa couplings, and fix the sparticle masses with the boundary conditions 
(0) and (|^). We next go down with the 6x6 mass matrices for the squarks and sleptons 
to find the sparticle spectrum. In running the RGEs, we include the one-loop correction 
to the running bottom quark mass. Ami,, which involves the contributions coming from 
gluino— bottom-squark loop diagram and chargino— top-squark loop diagram, and is given 



by EI 



Arrih = — Afofi/itan/? 



2a, 
3n 



MgI{ml,ml,Ml) + -^AJ {ml, fx 



(4) 



where the integral function I{a, b, c) is given by 

abln{a/b) + bcln{b/c) + caln{c/a) 



I{a, b, c) 



(a — b)(b — c)(c — a) 



(5) 
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and Mg is the gluino mass and (m^. ) is the bottom-squark (top-squark) eigenstate masses, 
respectively. Our convention for the sign of the Higgs boson mass parameter fi follows that 
in Ref. 0. 

Calculation of 6 — 57 amplitude involves the coefficients of short distance photonic 
and gluonic operators c^i^Mw) and csi^Mw)- The SM calculations of C7 and cs with QCD 
corrections to two- loop order are given in Ref. ||l9l. Calculations of next-to-leading order 
agree with the previous calculations while reducing the theoretical errors Various 
super symmetric contributions for 6 — > 37 are given in a generic form in Ref. |j21[]. The 
branching ratio for 6 — *• 57 is given by 



BR(6 ^ 57) 



TT I{mc/mb)[l - ^as{mb)f{mc/mb)] 



BR{b^ceu), 



(6) 



where r] = as{Mz) / cxsimb). I is the phase-space factor given by I{x) = 1 — 8a;^ + 8a;® — 
— 24:X^\n.x, and / is the QCD correction factor for the semileptonic decay given by 
f{mc/mb) = 2.41. and Ag are the coefficients of the effective 657 and bsg penguin 
operators evaluated at the scale Mz, respectively. BR(6 — ceu) denote the branching ratio 
of the semileptonic decay b — >■ cei'. 

The supersymmetric contributions to the muon anomalous magnetic moment are es- 
sentially coming from neutralino(x°)-smuon(/i) loop diagram and chargino(x"'')-sneutrino(z/) 
loop diagram as follows ||^,|^ : 



6a 



SUSY 



6 a,, + 6a]1. 



(7) 



Here 6a^ and 6a^ denote the contributions from the diagram and x^'^ diagram, 

respectively, and are given by 



6a, 



E 



A.a 



(iV^iV^ + iV^iV^)(l - Qx,^ + 3x1 + 2a:L - Q^a Inx. 



6m?^(l -Xi, 



:N^^N,'l{l-xl + 2xi^ Inxi^) 



6a) 



E 



m, 



I L 

3m-+ 



3ml{l-xi)^ 



{CtC[^ + Cl'Cl') (1 - hi - 3x^ + \x] + 3x, In 



Xi 



mr- 



where (i, / = 1, 2; a = 1 — 4) 



. 4 1 2 2, 



(8) 



(9) 



Xi. 



N. 



-^{U%^mu + V2gy{U'^] 



la 



Vl V2 v2 



m 

Vl 



V2' 



(10) 
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(72 and gy are the gauge couplings of SU(2)l and U(l)y, respectively, and vi is the vacuum 
expectation value of the Higgs boson Hi. m^, "m^o, m^+, mfx, and mc; are masses of the 
muon, neutralino, chargino, smuon, and sneutrino, respectively. and f/^^^ denote the 
neutralino and smuon mixing matrices, and U'-' and denote the chargino mixing matrices. 

We use our calculated mass spectrum and couplings to calculate the rate for 6 — > 57 and 
6a^^^ . For fixed values of tan/5, n and sgn(/i), both the branching ratio for b sj and 
6a^^^ as well as and the weak gaugino mass M2 are calculated, as the values of M and 
A vary. Then the bounds on the branching ratio for b ^ sj and (5a^^^"^ are translated into 
the bounds on values of M2 and in the — M2 plane for fixed values of tan/?, n and 
sgn(/i). From Eq. (Q) one can see that M2 is directly related to A, since g{x) ~ 1. Bounds 
on other sparticle masses can be easily deduced from a bound on M2, owing to the relations 
Eqs. (D and ©. 

In b —>■ sj decay, the contributions to the total decay amplitude are coming from the 
W loop diagram, charged Higgs boson loop diagram, neutralino loop diagram, and gluino 
loop diagram. It has been pointed out that the neutralino and gluino contributions to the 
amplitude are less than 1 % in the whole range of parameter space The charged Higgs 
boson loop contribution adds constructively to the W loop contribution, while the chargino 
loop contribution can be constructive or destructive to the W loop contribution, but is 
generally much smaller than the charged Higgs boson loop contribution. We use the new 
CLEO bound on the branching ratio for 6 — 57 decay in order to obtain constraints on the 
parameter space of the GMSB models. 

The bound on the supersymmetric contributions to is given by —71 x 10~^° < 5a^^^"^ < 
207 X 10~^° at 90 % CL. This bound is obtained by the difference between experimental 
value and theoretical prediction of a^. The new E821 experiment at Brookhaven is expected 
to improve the experimental determination of to the level of 4 x 10^^" Since the 

complete two-loop electroweak contribution in the SM to is = 15.1(0.4) x 10~^° 
and the supersymmetric contributions can be as large or even larger than the , it is 
expected that the new E821 experiment will be possible to test both the SM electroweak 
and supersymmetric contributions [Q-§|. 

In Figs. 1—8, we display the bounds obtained from the branching ratio for 6 — 57 and 
^qSusy ^Yie — M2 plane for either sign of /i, for tan/5 = 10 and 60, and for n = 1 and 
3, respectively. Solid lines represent the bounds from the branching ratio for b s'-f and 
dot-dashed lines describe the bounds from ^a^^^"^. Figures 1 and 2 show the bounds on M2 
and \n\ for tan/? = 10 and n = 1, and for positive and negative n, respectively. The region 
surrounded by the solid line is allowed by the CLEO bound, while the upper region of the 
dot-dashed line is allowed by the present bound on a^. In the case of Fig. 1, the constraint 
from b s'-f decay is clearly much stronger than that from a^. We find M2 > 248 GeV and 
fi > 626 GeV. Small values of M2 lead to unacceptably large contribution to the branching 
ratio for b ^ s'j, while large values of /z raise the problem of fine-tuning and are generally 
constrained by the lower bound on the stau mass. In Fig. 2 we see the constrains from both 
b ^ S'y and are complementary. By combining the bounds from the both, we can obtain 
much stronger bound on M2 and in particular, low values of which would have been 
allowed are excluded. We find M2 > 210 GeV and \fj,\ > 505 GeV. 

In large tan/? case, we find that the bound from either 6 — > S7 or is more stringent 
than that in small tan/? case, and most region in the — M2 plane is excluded. For 
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tan/3 = 60 and /i > (Figs. 3 and 6), the allowed regions from each of 6 — S7 and do 
not overlap, even though a possibility exists that they might overlap for unacceptably very 
large values of /x. Thus, this case is excluded, while it would be allowed if one considered 
only either 6 — 57 or as in Refs. |]6|,pJ]|. For tan/3 ^ 50 and /i > 0, the allowed regions 



from each of 6 ^ 57 and overlap allowing limited regions in the parameter space. For 
tan/5 = 60 and < 0, the supersymmetric one-loop correction to bottom quark mass leads 
to unacceptably large value of rrib. In other words, in order to have a correct value of m^, 
the negative sign of /i is not physically allowed for large tan/? in our analysis. It can be 
qualitatively understood from Eq. : for large tan (3, the correction Am^ is large, and for 
negative fi, Arrify gives positive contribution leading to an incorrect value of rrify. Thus, by 
inclusion of the correction Am^, we exclude the case of large tan f3 and ^ < 0. This result is 
different from those in the previous works : in Ref. [|lT| the one-loop correction Am^ was not 
taken into account and the case of large tan f3 and /i < was favored without any constraint 
from b —>■ 57. Also, in Ref. 0, the case of large tan f3 and /x < is allowed with a constraint 
on the value of M2. 

For n = 3, the constraints from each of 6 ^ S7 and are more stringent. But after 
combining the constraints together, we find that for tan/3 = 10 and /i < 0, the resulting 
limits are lower (i.e., M2 > 182 GeV and > 320 GeV) than those in the case of n = 1, 
while for tan/? = 10 and /i > 0, the resulting limits are still higher (i.e., M2 > 357 GeV and 
fi > 648 GeV). For n = 3, large tan/3 region is almost ruled out due to the same reason as 
the case oi n = 1. 

In Figs. 7 and 8 we plot the bounds on the sparticle masses, obtained by this combined 
analysis of 6 — 57 and function of tan /3 for positive fj, and for n = 1 and 3, 

respectively. The plots are displayed for up to tan /3 ~ 50, since the region corresponding to 
tan/3 ^ 50 is ruled out. The lower bounds on the sparticle masses increase monotonically 
as tan/3 does. For n = 3, the lower bound on each sparticle mass is higher than that for 
n = 1. (We show the plots in the positive fi case only, because for negative n the allowed 
parameter space is restricted to only relatively small tan/3 region (tan/3 < 20) in order to 
have a correct value of nib in our analysis.) 

Some comments concerning the non-leading order (NLO) corrections to b sj are in 
order. It has been pointed out by Kagan and Neubert ||2^ that NLO corrections lead to a 
10% error in predicted theoretical values. Ciuchini et al. find that, when the leading 
order (LO) corrections cancel, contributions from NLO corrections become significant; this 
happens when the masses of charginos and right-handed top-squark are very much higher 
than those of other squarks and the gluino. In our parametric space such a cancellation of 
LO corrections is unlikely. There could be other circumstances leading to the cancellation 
of LO corrections, and this is under investigation. 

In conclusion, we have investigated both the 6 — >■ S7 decay and the anomalous magnetic 
moment of muon together in MSSM with GMSB. We have used the new CLEO bound 
on the branching ratio for b —>■ sj and the present experimental limit on to constrain the 
parameter space of the GMSB models. We have presented bounds on supersymmetric parti- 
cle masses as a function of tan /3. This combined study has led to much stronger constraints 
on the parameter space of the model than those from either 6 — 57 or a^. In particular, 
the region of large tan /3, which would have been otherwise allowed, is ruled out or severely 
constrained, depending on the sign of fi. With the inclusion of the supersymmetric one-loop 
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correction to b quark mass, we have found that the region of large tan /3 and negative jj, is 
physically ruled out in order to give a correct value of m^. The present experimental data on 
the decay 6 — > 57 leads to more stringent constraints than those from in a broad region of 
the parameter space. The anticipated precision level of 4 x 10~^° in determination of in 
the Brookhaven E821 experiment would constrain the parameter space much more severely. 

We thank Jim Smith and B. Dutta for helpful discussions. This work was supported in 
part by the US Department of Energy Grant No. DE FG03-95ER40894. 
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FIGURE CAPTIONS 



Fig. 1 : Limits on the weak gaugino mass M2 vs |//| for tan/3 = 10, // > and n — 1. Units 
are in GeV. The soUd line represents the bound from the branching ratio for b ^ 
(the region surrounded by the sohd hne is allowed) and the dot-dashed line represents 
the lower bound from a^. Note that M2 and are calculated and not independent. 
The value of M2 increases as the value of increases in the allowed region. 

Fig. 2 : The same as Fig. 1, except < 0. 

Fig. 3 : Limits on the weak gaugino mass M2 vs fi for tan ,5 = 60, n > and n = 1. Units are 
in GeV. The solid line represents the bound from the branching ratio for b 57 (the 
region surrounded by the solid line is allowed) and the dot-dashed line represents the 
lower bound from a^. Note that M2 and are calculated and not independent. The 
value of M2 increases as the value of |//| increases in the allowed region. 

Fig. 4 : The same as Fig. 1, except n = 3. 

Fig. 5 : The same as Fig. 2, except n — 3. 

Fig. 6 : The same as Fig. 3, except n = 3. 

Fig. 7 : Bounds on the sparticle masses (in GeV) as a function of tan (3 for > and n = 1. 

The solid line represents the lower bound on the gluino mass, and the dotted and 
dot-dashed lines represent the lower bounds on the stop {rrif^ and m^^) and sbottom 
{mi^ and m^^) masses, respectively. 

Fig. 8 : The same as Fig. 7, except n = 3. 
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